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£h ■ Abstract 
£3 ■ 

The method for calculation of the correlation functions of the Ising- 
type systems with short-range interaction and with arbitrary value of 
spin is developed within cluster approximation. For the Ising model 
(spin S z = ±1) the expressions for pair and ternary correlation func- 
tions within two-particle approximation in q-space are obtained for 
the hypercubic Bravais lattices. In the ID case the exact expressions 
for them in the site space is obtained as well. On the basis of the 
Glauber equation within two-particle cluster approximation the longi- 
\ tudinal dynamical susceptibility x(cL E) is found. In the ID case and 

in the absence of external field the expression for x(q, E) is exact. For 
the Emery-Blume- Griffiths model (S z = —2, 0, 2) within two-particle 
approximation the pair correlation functions are calculated. The four- 
particle cluster approximation is used for calculation of static suscep- 
. tibility x(q) of KD2PO4 ferroelectrics. 



1 Introduction. 



c3 

a 

The wide class of materials is known which is described on the basis of the 
pseudospin Hamiltonian. In particular these are magnets with different val- 
ues of spin M [1-3] (in our notations S z = -M, -M + 2, . . . , M - 2, M), fer- 
roelectrics with hydrogen bonds [4,5] (M = 1), multicomponent alloys [6,7] 
(M + 1 corresponds to the number of components), lattice gas [8] (M + 1 
is the number of atom states per site). In the works [9-12] the method 
of investigation of pseudospin system based on the introducing the refer- 
ence system (RS) was developed. This approach makes possible to take 
into account adequately both the short-range and long-range interaction. 
RS contains besides the low-dimensional short-range {D = 1,2) part of the 
Hamiltonian the terms taking into account the other types of interactions 
in the molecular field approximation (MFA). The fluctuating corrections 
to MFA for these types of interactions are taken into account on the basis 
of expansion of the physical responses (free energy, correlation functions, 
temperature Green functions) in terms of loop diagrams or in terms of 
many-tails diagrams. These corrections contain the physical responses of 
RS. The calculation of the latters is the separate problem of the referense 
system approach. 

In the series of cases for the description of pseudospin system's prop- 
erties it is enough to take into account the long-range interaction in MFA. 
Now the reference system problem is to be considered separately. However, 
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only a small number of low-dimensional problems may be solved exactly. 
The exact solutions are obtained mainly in the case of ID quantum models 
or ID and 2D Ising-type models (ITM) (the Hamiltonian of ITM contains 
only S z component of the spin). Among approximate methods of investi- 
gation the cluster method is sufficiently effective [13-31]. It is necessary to 
distinguish the cluster expansion method (CEM) and the cluster variation 
method (CVM). In the case of CEM the cluster expansion is constructed 
for the free energy [13] whereas in the case of CVM one carries out the 
cluster expansion of entropy [26]. It can be shown that in the first order 
of cluster expansion both methods give equivalent results for all responces 
of pseudospin system. This approximation we shall call the cluster approx- 
imation (CA). Let us note that in the first works within the framework of 
CVM the results were obtained on the basis of combinatoric approach [25] 
whereas within CEM the results was obtained on the basis of many-particle 
generalization of MFA [14] . 

Unfortunately the cluster approach was used mainly for investigation 
of thermodynamical properties of pseudospin systems. We know only sev- 
eral papers [32-43] where the cluster approach was used for calculation of 
reference system's correlation functions (CF) (see [32-37]) and Green func- 
tions [38-43]. In the works [32,33] the equation for pair CF (Ornstein-Zernike 
type equation) of ITM was constructed artificially. The dipole-dipole inter- 
action which has the nonanalytical properties at q — ► was taken into 
account in MFA. The concrete calculations were carried out only in the 
case q — > 0. Let us note that only in paraphase the obtained results corre- 
spond to cluster approximation. In the work [34] within the frames of CA 
for ITM the method for calculation of pair CF for RS, when the long-range 
interaction is absent, was suggested. This method is applicable in the case 
T > T c . Here analytic and numerical investigation of ^-dependence for CF of 
some models was performed. The Ising model (IM) (IM corresponds to ITM 
with M = 1) on square lattice within the frame of two-particle cluster ap- 
proximation (TPCA) and four-particle cluster approximation (FPCA) was 
considered. The model of ferroelectrics KD^PO^ was studied too. In the 
paper [35] for IM within the frame of cluster expansion method the method 
which allows to find the correlation function of arbitrary order for RS within 
cluster approximation was suggested. For IM on hypercubic lattices within 
TPCA the pair and ternary CF were found. Here it was shown that for 
ID system these results coincide with the exact ones. In the ref. [36] it was 
proved that TPCA for IM yields exact results for all characteristics of the 
system in the case of tree-like lattice. Moreover in the case of square lattice 
within CEM the influence of the higher order terms of cluster expansion 
on thermodynamic characteristics was studied numerically. The pair CF 
in FPCA was also calculated. In the paper [37] CEM was generalized on 
the Ising type systems with arbitrary value of spin and detailed numerical 
investigation of the Blume-Emery-Griffits model was performed. 

In the works [38-41] the dynamics of ferroelectrics with hydrogen bonds 
with accounting of tunneling was investigated within the frames of clus- 
ter approach. Here the method of two-time Green's function was used. 
The equation for them was uncoupled in the spirit of Tyablikov approxi- 
mation with respect to long-range interaction. At such consideration the 
short-range intercluster correlations were not taken into account and intr- 
acluster Green's functions appear to be connected only through long-range 
interaction. This leads to uncorrespondence of expressions for static sus- 
ceptibility which follows (at q = and cj n = 0) from the expressions for 
dynamic susceptibility with expressions which are obtained within CA. In 
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the works [42,43] a method was developed which within the frames of cluster 
approach allowed to obtain selfconsistent results for thermodynamic char- 
acteristics and temperature Green's functions. Within TPCA the reference 
Green's functions were calculated for the Heisenberg model [42] and the 
Ising model in a transverse field [43]. 

In the present work on the basis of cluster approximation the pseu- 
dospin systems which are described by Ising type models, are investigated. 
In section 2 the statement of CEM for ITM with arbitrary many-particle 
short-range interaction is presented. In section 3 within TPCA the ther- 
modynamic, correlation functions and relaxational dynamics of IM (M=l) 
on hypercubic lattice are investigated. In section 4 in the frames of TPCA 
the Ising type model with pair interaction and arbitrary value of spin is 
considered. Section 5 is devoted to investigation within TPCA of the Ising 
type model describing the ferroelectric properties of crystal KD2PO4. 

2 Statement of the cluster expansion method. 

In this section we shall consider the cluster expansion method for the Ising- 
type model with arbitrary form of the Hamiltonian. In the following sections 
we shall apply present method to certain models. The Hamiltonian of IM 
with arbitrary values of spin can be written in the following form (Sf = 
Si = —M, —M + 2, ... ,M — 2, M) 

-f5 L H= L H({h^}) = H({h^}) + ±^2^^^), (2.1a) 
H({h { ) }) = E WSJ + E K a( S i, Sj) + W({S}), (2.1b) 

* (ij) 
M MM 

hi{Si) = E hV>S?; Kyi*, Sj) = E E 4 W ° W > ( 2 - lc ) 

H=l fi=l v=l 

M M 

Jij(SiS j )=Y,Y, j ir )s ? s j- ( 2 - id ) 

fl=l v = \ 

The one-site part of the Hamiltonian hi{Si) describes interaction of pseu- 
dospins with the fields of different type. The two-site part of the Hamil- 
tonian contains the short-range Kij(SiSj) and long-range Jij(SiSj) pair 
interactions. The three-site, four-site etc. short-range interactions are in- 
cluded in W({S}). The pseudospin system with the Hamiltonian H({h^'}) 
is called the reference system [9-12]. We shall consider the RS which con- 
tains only the interaction with the nearest neighbours. In particular, for 
pair interaction we can write 

K (^u) _ K {^,v)^ . _ _ / 1, i is the nearest neighbour of j, ft) , 
~ A nt i> n v ~ \ 0, in opposite case. ^- L > 

In the case of the Bravais lattice the indices i,j denote the sites of the 
lattice. The case of the structure with sublattices we shall consider in sec- 
tion 5. In the present paper the pair long-range interaction is taken into 
account in the molecular field approximation (MFA). Now the .F- function 
(the logarithm of partition function) of the system with the Hamiltonian 
(2.1) can be written in the following form 
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L H{h { ) }) = -\ E E 4T } L (^f ) L (^J) + HiC^}). (2.3) 

Here ^"({C^}) i s ^"-function of RS with the Hamiltonian containing the 
long-range molecular field A^: 

^({C«}) = lnZ({C«}); Z({C()}) = Sp w exp[ff({C«})]; (2.4a) 

C M = ^) + A (M). A (M) = ^ L W . (2 . 4b) 

Here the average L (S^) is taken for canonical ensemble and the density 
matrix is constructed with the Hamiltonian (2.1a). We shall calculate cor- 
relation functions (CFs) of RS (cumulant averages of products of operators 
) in the following way: 



fcf;{C ( - ) }) = W...sj 



* ^({C W ». (2-5) 



Similarly to (2.5) we shall obtain CFs L J r ^{^f l l ;{h^}) for system with 
long-range interaction. In this case the ^-"-function L J r ({h^}) is differenti- 
ated with respect to fields In the MFA these CFs can be connected 

with the CFs of RS. In particular for L (5f ) using the relations (2.3), (2.5) 
we obtain the equality 

L J™ (f ; {/»«}) = Hfl) = <Sf ) = (f ; • (2-6) 

If we shall differentiate (2.6) and take into account the functional depen- 
dence of C (m) on {S$ ), we shall obtain the relation for pair CFs. It has the 
following matrix form 

L fW{{h^}) = ^( 2 ) ({((')}) [i_ J^ 2 )({CW})] -1 • (2-7) 

Here the matrices A = are constructed with the help of power-type 

indices n, v and indices of sites i, j. In the homogeneous field h!f^ = 
the Fourier transformation diagonalizes the relation (2.7) with respect to 
the indices of the elementary cell. 

At first let us decompose lattice into clusters for construction of cluster 
expansion of F-function [36]. The figure 1 demonstrates decomposition of 
square lattice into clusters. 

Let us note that certain choice of lattice decomposition depends on the 
number of factors. In particular, this choice depends on structure and sym- 
metry of crystal, on form of interaction and cluster expansion order. In 
the present work we shall use only lattice decomposition into two-particle 
("link") and four-particle (tetrahedron) clusters. We shall denote the effec- 
tive field operator as Ripi(Si). It effects on site i from the side of cluster R, 
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Figure 1: The decomposition of square lattice into clusters ("link", 
"square") and demonstration of effects of cluster fields on site. Here "s" 
is the number of clusters containing present site, "b" is the number of clus- 
ters containing present link, "k" is the number of sites in cluster. 

which contains this site. The set of these clusters we shall mark 7Tj. Obvi- 
ously, Si fields effect on an arbitrary site i. After decomposition of lattice 
into clusters we pass from summation over sites to summation over clusters. 

]T £ RniSi) = EE WW); w({S}) = E w d{S}) , 

i R R feR R 

£*y(W = £ E ^-K ff ,{S f S r ). (2.8) 
(ij) R ff>eR u ff 

Here we shall consider the cases when Wr({S}) contains only spins from 
the cluster R. Taking into account (2.8) after identical transformation the 
Hamiltonian (2.1) is represented in the following way 

H = Y / H i (S t ) + Y / U R ({S}), (2.9) 

i R 

where 

H i (S i ) = h l (S l )+J2R^f(Sf), 

R 

Ur({S}) = - R<Pf(Sf) + (2-10) 
feR 

+ Y,r- K ff'( s f s f') + w d{S}). 

Then, using the form (2.9) and cluster expansion method (CEM) [13], the 
^"-function can be written in the following form 

T = J2 Fi + ln(exp (j2 Ur({S})) )o = E F i + E K i ■ ( 2 - n ) 

i \ R J i 1=1 
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Here we use the notations 



F t = lnZ t = \nSpe Hi ; Po = \{ Pl = J[ 



(2.12) 



Functions K\ contain sums over sets 1 = {R\ ■ ■ ■ Ri), which include noncoin- 
siding clusters Ri. 

K t = K(R 1 ,...,R l )=J2K(l), 

(Ri,-,Ri) {1} 



1 



1 



^(1)=E- E ^■■■—( U n\--- U nZ)o- (2-13) 

t/ l- I'm- 



ui=l v m =\ 



The cluster functions K(\), as it is well known, are expressed by means of 
L(m)-functions as well as functions L(l) are expressed by means of K(m) 
in the following way [13] 



l-i 



K(l) = ]T (-l) l - m £ L(m) = L{1) + ]T (~l) l - m E L ( m ) . ( 2 - 14a ) 

m=l mCl m=l mCl 

L(i) = E E * a) = *(m) + E E K ( m ) > ( 2 - 14b ) 

m=l mCl m=l mCl 

where L-function has the form (M(l)-moment function): 
L(l) = lnM(l) = F(l)- Yl F - 

- j 

E^n 



M(l) = (exp 



_n=l 



'0 • 



(2.15) 



Here we use the notations 1 for set of clusters, F(l) for JF-function, Z(\) for 
partition function and H(Y) for the Hamiltonian: 



F(l)=InZ(l) = 5p { s } exp{fr(l)}, 



H(l)= ]T H i+ll U Rn 



(2.16) 



n=l 



Arbitrary set of clusters 1 of present type forms certain graph (diagram) [36] 
on the lattice. We shall call graph to be unconnected one if it has parts which 
do not contain common sites. These parts will be independent statistically, 
when we shall carry out the averaging with distribution function pq. Let 
us separate arbitrary graph n into two graphs ni and 112 which have not 
comon clusters, but which may have common sites. 



ni ri2 n = niUri2, ni H n2 = M12 (2-17) 
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Here M12 is the set of sites, which is common for ni and ri2. Let us represent 
the sum £ mCn m the form X) m Em 2 i where mi C ni, m.2 C 112. Then we 
can write (after separation of items with mi = and rri2 = ) 

ni ni 

tf(n) = L(n)-L(m)-L(n 2 )- £ E E E K(ra 1+ m 2 ) . (2.18) 

mi=l m2=l miCni m 2 Cn 2 

The prime near the symbol of sum means that item at mi = n\, rri2 = ra 2 
is omitted. Let the following condition take place 

L(m 1 + m 2 ) = L(m 1 )+L(m 2 ), Vmi C ni , m 2 C n 2 . (2.19) 

Then we can carry out recurrent procedure for If (mi + m 2 ), which begins 
with mi = 1, m 2 = 1. For mi = 2i = li + l'i, m 2 = li we have 

if(li + l 2 ) = L(l 1 + 1 2 )-L(l 1 )-L(l 2 ) = 0, 

K(l' l + 1 2 ) = L(li + 1 2 )-L(li)-L(l 2 ) = 0, (2.20) 

K(2i + 1 2 ) = L(2i + 1 2 )-L(2i)-L(l 2 ) 

-if (li + 1 2 ) - if (l'i + 1 2 ) = 0. 

Similarly we can obtain 

if (mi + m 2 ) = 0, Vmi C ni , Vm 2 C n 2 . (2.21) 

The condition (2.19) is fulfiled for ni and n 2 being two unlinked parts of 
graph n (when Mi 2 = 0). The condition can be fulfiled also for some linked 
graphs, if we choose the fields R <pi by certain manner (see the following 
section or paper [36]). Let's note [36] that for set (2.14a) some L-functions, 
which correspond to unlinked diagrams are eliminated by L-functions of 
linked diagrams of certain type. Therefore there is a sum over linked graphs 
of certain class in (2.14a). 

In the present paper we shall use mainly the first approximation for 
L-function and the decomposition of the lattice into R-particle clusters 
(RPCA). As it follow from (2.11)-(2.16) the F-function in RCPA has the 
form [36] 

mc^}) = - E(^ - m({t { - ] }) + E fr({r< { - ] }) • (2-22) 

i R 

Here the notations for intracluster one-particle (i 7 ,) and R-particle (F R ) 
^"-functions and corresponding partition function Zj, Z R are introduced 

= m^ttcP}); = Spine**™-, (2.23a) 

Fr({r< {) }) = lnZfl^C }); Z r ({r< { - ] }) = Sp {s} e H ^ s \2.23b) 

The one-particle Hamiltonian Hi(Si) and R-particle Hamiltonian H R ({S}) 
have the form 

Hi(Si) = E M ? } Sf; ? } = Cl M) + £* 'r^; (2.24a) 
H R ({S}) = EfeRH j (S j ) + U R ({S}) = 

= E/ eR R(f(S f ) +E ff >eR vfpVfriSfSf.) + W R ({S}); (2.24b) 

rOW = flC^sjf ; rC { ; ] = Cf + Er>^r «<pf- (2.24c) 
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As we can see from (2.24c), the fields Rtpf 1 which effect on site / C R from 
cluster R are not contained in R-particle Hamiltonian. This is naturally due 
to the fact that all these interactions are taken into account exactly with 
the help of terms Vff and Wr. On the basis of Hamiltonians Hi, Hr we 
can introduce intracluster density matrices 

e Hi(Si) e H R ({S}) 
P m = — — ; p R ({S}) = — , (2.25) 

A ^R 

and intracluster CFs 

((S t yi(Sff 2 ...(Sf i yM)%, (2.26a) 

((S^ . . . {S%)"£ . . . (S fR f( R . . . {Sf E )^r pR . (2.26b) 

We shall calculate CFs (2.26) on the basis of F-functions (2.23) and rela- 
tions of (2.4)-type. Evidently, they do not equal zero only if the site indices 
(i or /) belong to present cluster (i or R). 

We shall find the system of equations for the fields Rip^ from condition 
of stationarity of ^"-function with respect to these fields. In the present 
paper we shall consider only RPCA. As it will be written in the following 
sections the system of NsiM equations in RCPA has the form 

W« = <sfU • ( 2 - 27 ) 

Here the average is denned with respect to density matrices (2.25). The 
system (2.27) is equivalent to the system of Nsi(M + 1) equations 

Pi (S i ) = Spn-iPR({S}), (2.28) 

which coincides with the exact relations between one particle and R-particle 
density matrices. Among Nsi(M + 1) equations (2.28) only NsiM ones 
are independent because Nsi conditions SpiSpn-i p r ({S}) = 1 are fulfiled 

identically. In the homogeneous field (hi = h, Q = £, R<p\^ = f^) (2.27) 
gives M equations for M unknown fields ip^ . 

3 Thermodynamic and dynamic properties of 
the Ising model. Two-particle cluster approxi- 
mation 

3.1 The main relations. Thermodynamic properties 

In the present section we shall consider the Ising- type model at M = 1, Si = 
±1 (Ising model- IM) in the case of hypercubic lattice. Now the Hamiltonian 
(2.1) contains only linear and bilinear forms of pseudospin operators and we 
can write it in the following way 

= E + E KijSiSf, K i3 = IK-,, . (3.1) 

* (ij) 

In TPCA the cluster contains only two sites (Ri r = (1, r)) and we can speak 
about effect on the site i from it's neighbour r. This fact leads to the change 
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of notations for the cluster fields: R<pi — > r <p±. Thus, taking into account 
(2.22)-(2.24), we obtain (sj = z, 6/// = 1) [35] 

^({C}) = (1 - ^)E^(Cl) + ^E^rUl ll Cr). (3-2) 
1 l,r 

The one-particle and two-particle cluster Hamiltonians have the form 

Hl(Si) = ClSi; Cl = Cl + Er'GTri r'¥>l > 

H 12 {S 1 ,S 2 ) = 2 Ci^i +i (2S2 + ^SiS 2 , (3.3a) 

2C1 = Cl + Er'gTri; rV2 r '^ 1 ' l( = 2 = C2 + Z)r'e7r 2 ;rVl r/</92- ( 3 -3b) 

Thus, the following expressions are obtained for the intracluster jF-functi- 
ons: 

Fi (Ci) = ^i (0) = ln{2 cosh Ci}; (3.4a) 
TP - T?W - 

= In {2e^[cosh( r Ci + iCr) + e~^ K cosh( r Ci - iCr)]} • (3.4b) 

If we take into account the long range interaction in MFA we must renor- 
malize the field Q: d = hi + A, = hi + J2j Jij(Sj)- Here hi is an external 
field, Aj a long-range-acting molecular field. 

In the present section we shall use the following notations 

ft.i R ({CD = (St, • • • S iR y = -A- . . . -A- ^({C}) , (3.5a) 

c£ } = 5 ij + E = Ac, 

&% = £r---w&= E r^!:2;i>3, (3.5b) 
(J-l) _ _<5_ 

4°(Ci) = ^|y [ ^ (0) (Ci) = (^ 1 , 

F 12 (2C1 li C2) - 9(2C> 9(i ^i2 -(^i^i)p 12 - (3.5c) 
Now taking into account (2.5) we easily obtain from (3.2) 

<Si> = A 1] ({(}) = -irHiQ) + E J^-^iC}) rrf} • (3-6) 
c*Ci r ^ 1 ' 

Here 
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— ^({C}) = (1 - z)F[ 1) + Y, F3 , (3.7a) 
-^-^({C}) = (1 - z)F§> + £ *iv " F$ . (3.7b) 

d r<Pl' 77 

Proceeding from stationarity condition for .F({C}) with respect to r ip\ 
(^-^({C}) = 0) on the basis of (3.6) and (3.7) we find the following 
relations 

mi = ( Sl ) = jf\{C}) = F^dCi}) , (3.8a) 
(Si) Pl = FW({Ci}) = Fg°>( r Ci |i Cr) = <5i> Wp • (3.8b) 

If the long-range interaction is absent (3.8b) gives the system of Nz equa- 
tions for Nz unknown variational fields r ipi({(}, K) while (3.8a) gives the 
expression for the average value (Si). The expression (3.8b) is equivalent to 
the relation between one-particle p\ and two-particle p± r density matrices: 

pi(Si) = Sp Sr pir(Si, S r ). (3.9) 

This equivalence can be proved if one takes into account that among 2Nz 
equations (3.9) only Nz equations are independent because NZ conditions 
Sps 1 s r Pir(S\, S r ) = Sps 1 pi{Si) = 1 are accomplished. It follows from (3.9) 
that 

*f(Ci) = {&)% = ((sj)^ = Fii\ci ii Cr) (3.io) 

The equation (3.9) can be written in the other form which is useful for 
diagram analysis 

M (lr) = Sp Sl s r (piPre Ulr ) = S PSr p r e Ulr . (3.11) 

Let us consider in expansion for ^"-function certain type of graphs which 
contains links that connect the site 1 with another part of the present 
graph [36]. In particular, this type of diagrams includes tree-graph, that 
is graph of links without loops. 



(3.12) 



Let us denote by v the set of all sites of graph n. Then for the moment 
function of graph (3.12) taking into account (3.11) we can write 

M(n + (lr)) = Sp {Sj} J] Pi e*p( £ U R )Spp r e u ^ = M(n)M(lr). (3.13) 
jev Ren 

Since M-function for the graph (3.12) is factorized, the corresponding L- 
function has the additive form (2.19). This fact means that for the graph 
of (3.12)-type the K function is equal to zero. It is shown in [36] that 
the differentiation jjr—^ of ^-function for tree type graphs also gives zero. 
But for i^-function of arbitrary graph this condition is not fulfilled (when 
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the index i in operator coincides with r in (3.12) and n does not 

contain links with free tails) . Thus if in the expansion of ^-function we shall 
consider only tree-graph and function r <p>i({Q},K) will satisfy the system 
(3.11) (or identical system (3.8b)) then the stationarity of ^-function with 
respect to fields will have been accomplished. In this case the expressions 
for /-"-function and (Si) coincide with expressions (3.2) and (3.8a) which are 
obtained in TPCA. It follows that for the lattice of the tree type, e.g. for 
ID chain, the TPCA gives an exact results for all quantities of the IM with 
Hamiltonian (3.1) [35,36]. 

Now we pay attention to thermodynamical properties of IM in TPCA 
taking into account the long-range interaction in MFA. In the case of ho- 
mogeneous field extracting (3 we can write 

r ifi = (3ip; r<Pi = z/3(p; & = PC = (3h + PJ m; (3.14a) 

r 

Ci = C = P(C + zip); rCi = P~C = P[C + (z- %]; (3.14b) 

where Jo = J(q — > 0). From (3.2), (2.23) and (3.14) for the free energy 
per site f(h,K) we obtain the following result 

p L f(h,K,J) = ij m 2 + /3/(C,K), 
-f3f((,K) = jffiC) = (1 - z)]nZi(Q + f lnZ 12 (C~). (3.15) 
Here the following notation is used 

Z 1 (C) = 2cosh(/3C); a = e~ 2 > 3K ; 
Zi 2 (Q = 2e^[cosh(2/^) + a] . (3.16) 
The system of equations (3.8) for m and ip has the form 

m = tanh[/3(C + zip)}; zip = -( + {2p)' 1 In 73^; (3.17a) 

tanh[/3(C + zip)] = sinh(2/3C)[cosh(2/3C) + a}' 1 . (3.17b) 

In the case Jo = the system (3.17) leads to single equation for <p(/3h, (5K) 
(or for m((3h, (3K)) if we exclude m (or p) from (3.17a). In the ID case 
(J = 0, z = 2) the equation (3.17b) for <p can be solved exactly, and 
for m and / from (3.15) and (3.17a) we get the generally known exact 
results [44,45] 

-Pf = K + lnfcosh PK + r], (3.18a) 



m = SinhW ; r = J a 2 + sinh 2 (P~h) . (3.18b) 
r v 

The equation (3.17b) can be used to obtain expressions for It 

(2) 

makes possible to find the pair correlation function ((SiSj) c )q = F q at 
q = and longitudinal isothermic susceptibility x zz ( c i. = 0) = \% z f° r LM: 
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Here d is the dipole electric (or magnetic) moment per one cell, v a is the 
volume of elementary cell, 

(J rf») r i = £ [a + cosh(2/3C)] 2 _ _z-l (3 _ 2Q) 
2 l + acosh(2/3C) l ~ m 
In the ID case from (3.20) using solution for (p one obtains the expression 

(^ 2 ))-! = r 3 [a 2 cosh^/i)]" 1 . (3.21) 

On the basis of (3.15) one can easy obtain other thermodynamic quantities 
of IM in TCPA, in particular, entropy and specific heat. In the case h = 
the equation (3.17) for ip (or m) always has the solution 92 = (m = 0) and 
at T < T c (below the phase transition temperature) these equations have 
also non-trivial solution p / (m / 0). As it was shown in [36], in TCPA 
for IM at K > the phase transition (PT) of the second order takes place. 

The equation for T c is found from divergence condition for L J 7 ^ 

z exp[-2/3 CJ FC] - z + 2 = 2(3 C J . (3.22) 

At Jo = from (3.22) one obtains well known result for T c of IM in 
TCPA [44,45] 

z 



k B T c = 2K 



i -1 



la- 



in the ID case (z = 2) from (3.23) we find that T c = 0. Thus for ID IM the 
PT is absent at temperatures above zero. 

3.2 The correlation functions 

In the present section we are going to consider briefly the suggested method 
for calculation of the correlation functions of IM in the cluster approach. 
We are going to find the explicit expressions for the pair and ternary CFs 
in the wave vector space. The results and the method of calculation of CFs 
in the site space in ID case are presented in the work [35] . 

On the basis of relations (2.5) and (3.8a) CFs can be expressed with the 

help of derivatives of fields Q ( se e (3.5b)) with respect to Q. We present the 
general expressions for pair, ternary and four-spin CFs in the wave vector 
space [35] 

^ 2 >(qi)=.F< 2 >C (1) (qi), (3.23a) 

^• (3) (qiq 2 ) = i ?(3) C {1) (qi)C {1) (q2) + i ?{2) C (2) (qiq 2 ), (3.23b) 
^ {4) (qiq 2 q3) = ^ {4) C (1) (qi)C (1) (q2)C (1) (q3) + 

+F (2) C ~(3) (qiq2q3) + F (3) [c( 1 )(q 1 )C( 2 )(q 2 q3) + 

+C {1) (q2)C (2) (qiq 3 ) + C (1) (q 3 )C (2) (qiq 2 )] , (3.23c) 

F (2) = i_ m 2. F (3) = _ 2m i?( 2 );F( 4 ) = -2F( 2 )(1 - 3m 2 ). (3.23d) 

Thus it is necessary to construct the equations for C^(q), C^(qi>q2) 
and C 3 Hqi>q2>q3) ( on the basis of (3.8b)) for the computation of ^ 2 ^(q), 
^(qi,q 2 ) and (q 1; q 2 , q 3 ). 
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Now we are going to consider the equation for (^(q). Differentiating 
(3.8b) with respect to C2 one obtains 

P (2)T(1) _ (20) (11) ~(1) 



^1 Cl2 - ^lr rCl2 + b \r lCr2 ) (3.24a) 

7(1) 
rS>i2 



?V) J 1 ). ?0)-?0) ,.(1) /"}9/lM 

"S12 - (.12 ~ r^l2 ' !Sr2 ~ V2 ~ W r 2 ■ (3.24DJ 



Taking into account (3.24b) and (3.10) the equation (3.24a) and the same 
equation with replaced indices l<->r constitute the system of equations for 

the fields an d 1^2 '■ This system can be written in the form 

u 'r)(J)-tefh-fP (a25) 

From (3.25) one can easy express r ip^ in terms of Taking into con- 
sideration the relationship (3.5b) between and J2r rf^j one can obtain 
the equation for (ff. 

SiiCS =fe + E <*Trl CS } • (3-26) 
r 

Here the following notations are introduced 

©ii = 52 " — j — )f lr f rl . _ 1 _ f lr f rl = d rlr . (3.27) 

Let us note that (3.26) is the equation for "correlation function" (f^ of 

(2) 

the Ornstein-Zernike type. If we take into account the relation T\j = 

F^Cff we obtain the equation of the present type for pair CF . In 

the works [32,33] the equation of this type was constructed artificially. This 
fact leads to unselfconsistent results for static susceptibility obtained from 
CFs and from free energy in the case of ordered phase. For uniform field 
the following relations take place 

fir = /7ri r ; din = d = 1 - f 2 ; On = & = ^^"^' i (3.28a) 
?2Q _ l + a 2 + 2acosh2/?C_ pU _ 1 - a 2 



F zv = 11 S ' F 



[a + cosh 2/?C] 2 [a + cosh 2f3(] 2 

/ = — ^ - 1 ^7§j = ' (3.28b) 

1 + a 2 + a cosh 2/3C cosh(/?C) + r 

Performing Fourier-transformation for (3.26) and taking into account 
(3.23a) we obtain the expression for the pair CF of IM in TPCA 

tV) = F m i±I = (°) (o on) 

fo) 1 - (* -1)/ +(1-/)- 1 /z0(q) l + <DzG(q)- l " ' 
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Here the following notations are introduced 



(2) F^(l + f) _ fjl-f)- 1 , 

-l-(z-l)/' *-!-(*-!)/' (3 - 30a) 

1 2 D 



9(q) = -[vr(O) -Tr(q)] = - £ sin%/2); 
7r(q) =^vr(r)e- iqr = 2^cos % . (3.30b) 



i=i 



Let us note that equation (3.26) in site space is exact for the tree-type lat- 
tice. But for Fourier-transformation (3.30) we use the translational symme- 
try. Thus the expression (3.29) is exact only for ID case. It leads to the well 
known results, which we write in both wave and site spaces (q = qi) [46-48] 



( 2 )_ COsh PC a 2 r 1 _ /oov/fi^-n , , 

•T(q) - 9 i (1 2^ ■ 2 /o ' WO^n ~ \ r ) J ■ [6.61) 

r z + (1 — a z ) sin g/2 



The expression for (SoS n ) is also obtained in the work [35] in another way. 

(2) 

Now we are going to consider how to obtain the equation for Q^j . It will 

give us the possibility to calculate the ternary CF. Differentiating (3.26) one 
finds 

M) tit 

ft 7( 2 ) _ hrtJrl + HrJriz 

U HS123 — ~ / , ^2 ^12 + 

r a lrl 

+ E ^V^ 13 ^ + E T^CS • (3-32) 

r °lrl r Ulrl 

Here the following notations are used 

fill = 4rflr, fill = ^frl • (3-33) 

The derivatives (3.33) one can easily compute: 

f W _ f (oi) .-i r ;(i) - ;(i) 

Jlr-3 — /lr a lrllV3 JrlS: 



13 J' 



where 



ffi = «[^-/irCS J ]- (3-34) 



f (01) _ 9 f _ _ 9f ff(01)_ f p(10)i 

/2 0) = ^r/ir = 0; / r ( i 0) = -t/ rl = 0. (3.35) 

<9rCl <?lCr 
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Substituting (3.35) into (3.32) one can find the equation for ternary "CF" 

a(2) 

enCl^E-^CS (3.36) 

+E KOT + b^iOT + OT] - c lr mi ] } ■ 

r 

Here the notations are introduced 

A - [f( 01 ) tl f (oi) f 1,-3 
-K-lr — [Ji r j r l J r i Jlr\ a i r i > 

Blr = [fTfl ~ ffr^frlKA , (3.37) 
C\ r = [frl ^ fir ~ fir ^Hlrl ■ 

Taking into consideration (3.28) the equation (3.36) for uniform field can 
be written in the form 

©Ci ( 23 = ^E^CS (3-38) 

r 

where 

2m/ 2 ^_ 2m/(l + / + / 2 ) 

(TTTFrf' (i + /) 2 ^ • ( } 

Note that (3.38) is the three-particle analogy of Ornstein-Zernike equation. 
Performing Fourier-transformation in (3.38) and using the relation 

9 - 4(qj + q 2 ) = 1 (3.40) 

« C (1) (qi + q2) 

we get the following result for C^(qi,q 2 ) 

C (2) (qiq2) = {[z + 7r(qi) + 7r(q 2 )]A - C7r(qi + q 2 )} x 

xC (1) (qi)C (1) (q 2 )C (1) (qi + q 2 ). (3.41) 

Substituting (3.41) in (3.23b) for ternary CFF( 3 )(qi, q 2 ) we finally find [35] 

•^ (2) (qiq 2 ) = {^ (3) +F( 2 )[(z + 7r( qi ) +7r(q 2 ))A- 

-CV(qi + q 2 )]C (1) (qi + q 2 )} C (1) (qi)C (1) (q 2 ). (3.42) 
It is easy to see (invoking (3.39)) that at T > T c the relation ^(^(qi, q 2 ) = 

(2) (3) 

fulfils. For computation of Ci 2 3j ^123 ^ ^ s necessary to transform (3.41), 
(3.42) to site representation. The expressions for them (see [35]) will contain 

the sums over sites. In the ID case these sums due to the form (j^ = /' 1_2 ' 
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lead to the sums of geometrical progression. After cumbersome transforma- 
tion the following expression for ternary CF can be obtained [35] 

F^(j,k) = (SoSjStf = F^f k ; (W(j,k) = 2mf k (f - 1); < j < k. 

(3.43) 

It is important to note that J-^ (j, k) does not depend on intermediate index 
j (0 < j < k) . Let us note that the expression for J 7 ^ (j, k) which is ob- 
tained by us in the frame of TPCA coincides with the result of works [47,48]. 
In these works the explicit expressions for mentioned above quantities are 
written only at k = j + 1, j + 2. In the work [35] using the solution for the 

fields r <p( 2 \j,k) the expression (3.43) is obtained immediately in the site 
representation. 

3.3 Relaxational dynamics. 

In the present subsection we shall consider the relaxational dynamics of the 
IM. Our main task is the calculation of the longitudinal dynamic suscepti- 
bility of the system under consideration 

/+oo 
die-™ ^e^X^i), 

Xij(t,t') = j^\sc jt ,=o; mt = {Si) P (t) = m i + Sm it . (3.44) 

Here we consider the dynamical response of the system with respect to small 
time-dependent mechanical perturbation of the initial Hamiltonian in the 
form J2j SQ(t)Sj . For the system described by the Hamiltonian (3.1) we 

can calculate only the static susceptibility x(q, 0). The evaluation of dy- 
namic susceptibility can be accomplished if we shall take aditional to (3.1) 
items which contain other component of pseudospins or if we shall take 
into account the interaction of the IM with dissipative subsystem (thermo- 
stat) [48]. Further we shall investigate the system which contains becides 
the (3.1) with the substitution Q — ► Q jt = d + ^d,t also the Hamiltonian 
of thermostat and interaction of thermostat with the Ising subsystem in 
the form ^2j a UjSj (a = x,y,z). Here L7" contains only the variables of 
thermostat. Within the framework of nonequilibrium statistical operator 
(NSO) method one can obtain [48-51] the kinetic equation (KE) for unary 
CF (Si) p uy Using some approximations [48] this KE can be lead to the 
Glauber-type form [52,46] 

[1 + nnu](S t ) p(t) = (F^(eu)) p( t) , (3.45a) 
d d 2 

^it = ^ + &/ Ti ~fa2> eit = Ci* + XI K irS r ■ (3.45b) 

Here Tj = 7" 1 is the bare time of relaxation (7, is the bare damping) which 
depends on the parameters of thermostat and on parameters of interaction of 
pseudospin subsystem with the thermostat. In addition to [48] we introduce 

also the bare resonance frequency uji (Q = uj~ 2 ). This is entirely in the spirit 
of phenomenological equation for damped oscilator (see [5]) with frequency 
Ui and with the damping 7^. The equation (3.45a) in the case = can be 
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obtained also if we shall proceed from the Glauber master equation for the 
density matrix of the system [52,46]. If we shall develop F^fat) as a set 
in S r one can obtain expression for (F^fat)) which contains distribution 
functions starting from unary (Si) p M up to the function of the power z 
) P (t)- The site indices of mentioned above functions belong to 
the certain coordinate sphere tti of the site i. In the following for average 
{F^(eit)) we shall restrict ourselves only by one-particle approximation 

(FU(e it )) P(t) *FU(( it ), (3.46a) 
Cu = Cit + A lt ; Ait = E r¥>« = E K ir v {Sr)t (3.46b) 

and two-particle one 

(& 1 He u )) Ht) n{Fl%& + K ij S j )) Plt) = 

= LiifCit, Kij) + P(jCit, K^) k (Sj) t , (3.47a) 

jCit = Cu - j<Pit = Cit+ E rfit ■ (3.47b) 

r(zTTi;r^j 

Here the following notations are introduced 



L(jCu, K^) = Lij j = \ E (jCit + K tJ a) , (3.48a) 

Z <r=±l 

P(jCit, Kij) = Pij, t = \ E (jCit + Kijo). (3.48b) 



2 .= ± i 



The field r ipn is taken to be the same in both approximations. The 
equations for r <pn (or for (Si)t) will be written below. The direct evaluation 
gives the following results (a^- = exp(— 2(3Kij))\ 



L sinh[2/?jCa] 
tJ,t cosh[2/3 jCit] + cosh 2(3 

_ 2ajj sinh[2/3 jCit] 

1 + a?- + 2aij cosh [2(3 j Cit] ' 

p = smh2(3K 12 
lht ~ cosh[2(3 jCit] + cosh2(3Kij 

1 + a?- + 2ay cosh [2/? j Cit] ' 



(3.49a) 



(3.49b) 



P — F^HjCit, iCjt) I f /q A Q \ 

" mwoc.. .(,.) k - c A (3 9c) 

In the last relation we take into account the expressions for intracluster 
CF F^ l \ F^ 0) (see (3.28a)). Substituting (3.46a) and (3.47a) at i = 1,2 
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in (3.45a) and writing two-site systems of equations in the matrix form we 
obtain i^m it = 1 (S i ) t ; k m it = k (Si) t ) 

TiQ.it 1 m i t = F^tCit) - 1 m it ; 

d t k m t = L t -V t k m t ; d t = {<%rAt}; (3.50a) 
By the direct computation one can obtain the following relations 

Fl° = Vf 1 L t ; A = 1 - Pi2,t P2i,t ; 

I F±2,t \ _ _J_ ( Li2,t + Pl2,t Al,t 

V F 2i,t J A V L ^ + P21 >* Ll2 >* 

Here the function F^t = F( 10 )( 2C.1t, 1(2*) already was used by us in the 
subsection 3.2. The system (3.50) can be represented in the form 

nnimu = F$ - 1 m ht , (3.52a) 
A -1 [ri fla fc m it + Tj % fc m jt ] = F$ t - k m i;t . (3.52b) 

The dynamic version of TPCA is constructed similarly to the static case. 
Thus the equality 1 mu = fc m^ gives the equation for the fields jipu- 

For the static fields Cu = Ci the left side of (3.52) is equal to zero and 
we obtain the system of equations (3.8) for the nonfluctuating part wij. Let 
us represent the effective fields in the form of sum of static and fluctuating 
parts. 

Cu = d + Kiu j<Pi = j<Pi + 5 jipu ; (3.53a) 
Cu = Ci + Ku\ jCu = jCi + 5 jCu ; (3.53b) 
5 jCu = Ku - 5 jtpu; ^2 S j(fit = 5(it - Kit ■ (3.53c) 

Let us substitute (3.53) in (3.52) and obtain two systems of equations. 

The first system of equations for rrii coincides with (3.8) and the second one 
for Srriit has the form 

(2) ~ 

TiQitSmu = F> ' dQt - Sm it , (3.54a) 
A7 1 i T i Ma Smu + Pij Tj Cljt 5m j t ] = 

= Fg 0) [SCit ~ 5 mt ] + 4 U) [SCjt ~ 5 np jt ] - 5m ht . (3.54b) 

Here the fact that in TCPA Smu = S 1 rriit = S k rriit is taken into account. 
Using (3.54a) the terms of type TiOJuSmu can be excluded from (3.54b). 
As the result we obtain the relationship which connects Sjtpu, Sufjt, Kit, 
Srriit. Just one more analogous equation can be obtained by permutation of 
indices i <-> j. From this system we express the field Sjipu via SCit and Sma- 



(3.51) 
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Using relation (3.53b) we obtain the equation which one can write together 
with (3.54a) in the case of uniform field in the following form (r^ = tg) 

F^SQ t = [1 + T G n t ]5m lt , (3.55a) 
1 + (z - 1)/ 2 F( 2 °)] SCu = 
= ^ (20) (1 - f 2 )$Qt + /E ^rSnirt . (3.55b) 

Substituting F^SCu from (3. 55a) in (3.55b) we obtain the closed set 
of equations for bran. Differentiating the last equation with respect to 
and using (3.44) we find the expression for dynamic susceptibility of the 
reference system 

*(q^) = i+or q | < v nH = Mi + — H. (3-56) 

1 + \1[u)t{c\) tg 
where x(q) = -^ (2) (q) (see (3.29), (3.30)) and 

r(0) 1 + (z - l)/ 2 

T(q) = TT^e^y r(0) = TG (i -/)[i- (,-i)/] • (3 - 57) 

In the ID case (z = 2, J(q) = 0) when the fields C are absent (C = 0) 
from (3.56) (at ( = 0) one obtains the expression which coincides with the 
well-known exact solution of kinetic equation (3.45) (see [46,53]). This one 
gives (a = exp(-2/3if), q = q{) 

a tg 1 + a 2 

x(q) = — 2 ^ ; T-(q) = 



a? + (1 - a 2 ) sin 2 g /2 ' 2 a 2 + (1 - a 2 ) sin 2 g/2 ' 

(3.58) 

In the work [54] the cluster solution of (3.45) {Q = 0) for x(q, u) at arbitrary 
value of z was presented only at ( = 0, q = 0, and this result coincides 
with (3.56). Let us now take into account the long-range interaction within 
MFA. Now, as before, the equation for ran = {Si) p (t) nas the form (3.45a) 
but we use the system (3.55) in which 5(u = Shu + J2k JikSmkt- We find the 
susceptibility k x{q,u) with the aid of derivative L Xij(t — t') = Smu/Shni at 
5hu = 0. This one leads to the relations of (2.7) and (3.56) type 

l i \ x(q,^) L x(q) (0 r n \ 

x q ' w = i TTV~( V = i ~T" c) ( \ l t \ ■ 3 - 59 

l-^(q)x(q,w) l + 0(w)^r(q) 

In the expressions for L x(q) an d L/ r(q) it is convenient to separate the part 
that is independent on wave vector 

Vq) = ■ Vq)= Lr(0) ■ 

*^ l+i$(q)' TW l+i$(q)' 

L cD(q) = L d»z9(q) + L x(0)[Jo - J(q)]; Jo = J(q - 0) , (3.60) 



where 



■ x(0) .^!!(Lt£ ;Mo .i±(^iM! ; , # _ ' 



2? ' v (l-/)2? ' (!-/)£' 

X> = l-(z-l)/- J F 20 (1 + /). (3.61) 
If q — > the expressions (3.60) coincide with the results of the works [55,56]. 
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4 Application of the two-particle cluster approx- 
imation to the Ising model with the arbitrary 
value of spin. 

In the present section we shall consider the Ising-type model with the arbi- 
trary value of spin M taking into account only pair interaction (W^dS*}) = 
0). The series of concrete results we shall obtain at M = 2 for the Blume- 
Emery-Griffits (BEG) model: 

H({C} = £[CiS + (>S?} + ^[KSiSj + K'SfS]} (4.1) 

i (ij) 

which was suggested in the work [57] for the investigations of phase tran- 
sition in mixture of He 3 -He 4 . This model corresponds to the system with 
the Hamiltonian (2.1) at M = 2 in the case K^ v) = S^K^n^, j\f v) = 
Sfu,J^. In what follows below we do not write the index 1 and instead of 
index 2 we use the prime symbol. Let us note, that at (' = —^zK' the 
Hamiltonian (4.1) describes the system with the bilinear interaction KSiSj, 
with the quadrupolar interaction K'(Sf — 8/3) (Sj — 8/3) and with the ex- 
ternal field C [2]- Let us note that thermodynamic properties of the BEG 
model were investigated by various methods [2] including TPCA. However 
the calculation of correlattors have not been carried out untill present time. 
We will computate here within TPCA the pair CFs for the IM with arbitrary 
value of spin [37] . 

As it was noted in TCPA we can use the relation (2.22)-(2.26) with 
Si = z, bff> = 1, W({R}) = 0, R = (12). Now the .F-function of the system 
under consideration takes the form 

= (1 - z) £ ^({8°}) + \ £ F lr ({ r C { 1 ) }, • ( 4 - 2 ) 

1 l,r 

The one-particle Hamiltonian and two-particle one in the case of BEG model 
read [37]: 

H 1 (S 1 ) = ( 1 S 1 + ('S 2 1 , 

Ci = Ci + J2 *<Pi; ([ = ([+ E *<Pi , (4.3a) 

r'67ri r'G7ri 

Hn(SiS2) =2 CiSi +2 (iSf +i (2S2 +i (2S2 + KS1S1 + K'SjS 2 , 
2C1 = Ci -2 2C1 = Ci -2 <Pi ■ (4.3b) 

Let us write the expression for intracluster jF-function of BEG model 

Fi(Ci.Ci) = ln2[e 4 ^ cosh2Ci + 1/2] , (4.4a) 
FuUi, 2C1 I 1C2, iC 2 ) = ln2{e^^ eK ' 
[e 4K cosh 2( 2C1 + 1C2) + e~ 4K cosh 2( 2C1 - 1C2)] + 
+e 4 < cosh2 2 Ci + e 4 < cosh2iC 2 + 1/2} . (4.4b) 
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Similar to (3.5) we shall use the notations for derivatives 



Ml jU 2 

U *2 



C ~(>i M = *gMi) (fii m) = s _i^_ (45) 

'12 '12 

Fl ({Cl }) -W¥»W¥ {{c }) ' 

F (fc M ^H'^+'« e ) ({r g-) } | = (4.5b) 

-^({,Ci () }|{iC^})- 



9[ r c?¥" 5[ r cT¥" aiic^] 1 * »[i3 e) ]'' 

Using the relations (2.5) and stationarity condition d , u] FilQ^}) = one 

can find the system of N{z + 1)M equations in N Z unknowns (<5f ) and 
NzM unknowns rV'f({C ( ' ) }, {^ ( ' } }) 

jrd) (f ;{C (D } ) = (Sf 1 ), = (S^ = ^({g }), (4.6a) 

^ i) ({Ci ( - ) }) = (^ i ) P1 = 

= (sr i w= J p 1 ( r |o) (Ui (0 }i{ic»})- (4.6b) 

When the long-range interaction is absent (4.6a) gives the system of equa- 
tions for r^i^ while (4.6a) gives the expression for unary correlation func- 
tion (5f) p . 

The system (4.6b) is identified to Nz independent relations between trial 
one-particle pi(S\) and two-particle pi r (Si, S r ) density matrices (see (3.9)). 
These relations and some other ones in which the quantities p\ , p\ r and U± r 
have not explicit form remain the same form as in the case M = 1 (see 
(3.11)-(3.13)). In particular, in the expansion of ^"-function the .fT-function 
for the graph of (3.12) type give zero contribution. Similar to the case 
M = 1 the TCPA gives exact results for responces of the IM with arbitrary 
M on the tree-lattice. Let us note that below we shall use the relation of 
(3.10) type which we can write in the form 

<(Si)* • . . {Sf) K ™Y P1 = <(Si)* . . . (S M ) K -% lr . (4.7) 

Further we shall present the explicit form for series of responses of BEG 
model in TCPA over short-range interaction and in MFA over long-range 
interaction [37]. For the free energy function per site from (2.3) and (4.2)- 
(4.4) we find 

P L f = § [W + J' q 2 ] + 0f(C, C, K, K') , (4.8a) 
-Pf=(l-z) In Zx(C, O + | In Z 12 (l () . (4.8b) 
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Here the notations are used 

Z 1 (C,C') = 2e 4 ^' cosh2/3C + l, 



Z 12 = 2e^ +^ K '+ 4K ) 

+4e 4/3C ' cosh 2/3C+1, 

( = ( + z<p, ~C = C + (z - l)<p; 

c = c+ V, C = C'+(z- w 



cosh4/?C + e- 8/3 ^ 



+ 



- h + Jom , 
C = h' + J'q 



m 



(Si), q = (St), Jo = J(q - 0), 4 = j'(q 



0). 



(4.9a) 

(4.9b) 
(4.9c) 



The system of four equations (4.6) for m, q, ip, ip' on the basis (4.4)-(4.5)can 
be obtained in the form 



m 



4e 4 ^' sinh2/3C 
Zi(CX') ' 
_ 8e 4 ^~' cosh 2/3C 

9 " Zi(C,C) ' 

e 4 ^'sinh2/3C _ 

^i(C, CO 

_ e m+^K'+^K) sinh4/^ + e^' sinh2/^ 

^12(C,C) 

e 4 ^"' cosh2/?C _ 

Zi(CX') 

_ e^( 8 ^ +16 ^+ 4y )(cosh4/3C + e~^ K ) + e 4 ^"' cosh2/?C 

^12(C,C) 



(4.10a) 
(4.10b) 



(4.10c) 



(4.10d) 



Let us note that numerical investigation of system (4.10) at different values 
of parameters of the Hamiltonian (4.1) in the case Jo = Jo = is carried out 
in the work [37]. Here it is also investigated the thermodynamic potential 
and static susceptibility of BEG model. 

The method of finding CF of IM with arbitrary value M is similar to one 
described in subsection 3.2. But now the matrix relation over indices fx, v 
takes place. One can obtain from the expression (4.6a) the matrix relations 

of (3.23) type connecting CF with derivatives from (jf^ . In the present work 
we restrict ourselves to finding only pair CF. Thus differentiating (4.6b) with 

respect to C^ 2 ^ we obtain the relation of (3.25) type 

rV?12 + fir l<f>r2 = flrCr2 > 
frl rV?12 + 1<£V2 = /rlCl2 ■ 

Here the following notations and relations are used 



(2) 



-1 



13 



(11). 



F 



(2) 



£.(20) 
ij 



(4.11a) 
(4.11b) 

(4.12a) 
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jtfO = {F^acT})};^ = {F^\{,ch I {4 } })} , 



c. 



H v 



(4.12b) 



So expressing the matrices r ipij with the aid of Qj and forming the closed 

equation for one can find at the case of uniform field the expression for 

C(q). On the basis of matrix relation of (3.23a) type we find the expression 
for pair CF ^ 2 )(q). It can be written in two forms [37] 



(i-/r 1 /^(qf 1 (i+/), 



(q) = (l-z) 



+ z 



p(2) + £(11) 



+ 



+ 



(4.13a) 



(4.13b) 



In the case of BEG model we shall present the explicit form of ma- 
trix elements (intracluster pair CF), which are included in (4.13) (F( 1+1 '), 
F( 1+0 ! 0+1 ))-nondiagonal elements and I (Z')denotes the number of deriva- 
tives with respect to unprimed (primed) field [37]: 



F (2+0') = ge 4/3C' J 2e W + cosh2/3( ?] Z 

(0+2 , } _ 32e 4 ^' cosh2/?C 

zt 

P (i+i') - P (i'+i) - 32e 4 ^'sinh2/3C , 

z\ 



l > 



(4.14a) 



F (l+0'|l+0') =8e 8/3C J e 16^' 



4e/ 3(8C+16X') sinh8/3K+ 



+8e 4/ ^ sinh Af3K cosh 8/3( + 
+e 4 ^cosh4/3C~-e- 4 ^ 



2 sinh 2 2p(\Z u 2 , 



(4.14b) 



F (o+i'|o+i') = a±_ e m ^(^+4^) [cosh4/3 ^-8^ ] _ 2cosh 2 2/3C ~l 
F (i+o'|o+i') = F (o+i'|i+o') = ™_ e w Le^' +WK 'hmh2[3(x 

Z 12 I 

x sinh 4(3K +e P(™K'+4K) s i n h 4/3C - 2 sinh 2/?C cosh 2/?C J • 

Let us note that the form (4.13) is the typical one for TPCA. Formula 
(4.13) gives expression for averaging over configuration of pair spin-spin 
CF [58] in the case of Ising mixture with the quenched disorder of replace- 
ment. Now the indices fj,, v number sorts of atoms; and intracluster CF 
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contain averaging over sort configuration. In the case of the one-sort quan- 
tum system (the Hamiltonian of such system contains different components 
of spin operators S a ) the form (4.13) gives expressions for the pair cu- 
mulant temperature Green functions in the impulse-frequency representa- 
tion [42,43]. Now the indices fx, v number the x, y, z component of spin 

operators; — > F^ 2 \uj n ), — > F^ n \uj n ) are expressed with the aid 

of intracluster Green functions. 

5 Thermodynamics and pair correlation functions 
of KD2PO4 ferroelectrics in four-particle cluster 
approximation. 

5.1 Statement of cluster expansion method. Intracluster 
correlations function. 

We shall consider a system of deuterons moving on 0-D..0 bonds in com- 
pound of KD<iPO{ type. This compound is related to orthorombic symme- 
try with the four molecules per West's elementary cell which is characterized 
by basic vector: |ei| = |ei| = 7.469^4, |e3| = 6.975A The unit cell of the 
Bravais lattice is composed of two neighbouring PO4 tetrahedra together 
with four hydrogen (a =1, 2, 3, 4) bond relating to one of them (A - type 
tetrahedron). Hydrogen bonds going to another (B - type) tetrahedron be- 
long to four nearest structural elements surrounding it. Their coordinates 
are ni, ni +r2, ni +r3, ni +r4. The projection of two neighbouring tetrahe- 
dra on the plane (ei xe2) as well as the effect of cluster fields on deuterons 
are presented in the figure 2. 

The Hamiltonian of deuteron subsystem taking into account short-range 
interaction and taking into consideration long-range interaction J aa (n — n') 
in MFA can be written as follows [32,33] 

H{{Q) = E (naSna ~ £ V R ({s}) , (5.1a) 
n,a R 

Vr(W) = V(n l ,n 2 ,n 3 ,n A ) = £ V aa'Sn a aSn a ,a' + 

a,a' 

+ (5.1b) 

lb 

~^~^rii+r2,n2^rii+r3,n^^ni+r4,n4 \ > 

Cna = h na + X na ', \ na = Jaa'( n ~ n )(S n 'a') ■ (5.1c) 

n'a' 

Here S na = ±1 describes the position of deuteron in double well po- 
tential which situates on ath hydrogen bond in nth cell of Bravais lattice, 
n a = n + r a where r a is the radius-vector of the hydrogen bond a in the 
present cell n. In (5.1b) the first term (the first product of Kroneker sym- 
bols) describes the short-range configurational interactions of deuterons near 
"A" tetrahedra, the second term describes the same interactions near "B" 
tetrahedra. Let us note that parameter J aa (n — n!) includes also an indirect 
deuteron interaction via lattice vibration. Further we shall use the following 
well-known notations 

y 12 = V 23 = V34 = V41 = ~Wi = V s , 
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Figure 2: Two neighbouring tetrahedra (A and B type) and lattice con- 
structed from four-particle clusters in projection on the plane (ei x e2). The 
tops of clusters coincide with the centres of hydrogen bonds. The effect of 
cluster fields on sites of cluster lattice is denoted by pointers. The sites of 
cluster Rq = R(R\) are numbered by indices io = 
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Via = V 24 = -e + ^W! = U s , 

$ = 4e - 8W + 2Wi; e = e a -e s , (5.2) 
W = ei-e s , Wi = e -e s . 

Here e s , ei, eo, e a are four different configurational energies for sixteen con- 
figurations of deuterons near [PO A f~ tetrahedra (see [4,5,32]). It is conve- 
nient to use the general index % = (n^ttj). Let us decompose pseudospin 
lattice into clusters. The cluster we shall take in the form of tetrahedron. 
Its tops coincide with the centres of hydrogen bonds. The projections of 
the same clusters on the plane axb is noted by shade line (see figure 2). Let 
us note that the lattice which is formed from the same clusters is similar 
to deformated square lattice with chess decomposition of this lattice into 
square clusters. Evidently that the structure of the Hamiltonian (5.1) agrees 
with that of (2.1). So we shall use the relations of section 2 (starting from 
(2.8)) at R = (1,2,3,4) when we shall perform the cluster expansion. In 
the present work we limit ourselves to four-particle cluster approximation 
(FPCA). Thus one may use relations (2.22)-(2.26). For reference ^"-function 
we write 

H{0) = - E m) + E F *(Ut}) . (5-3) 

i R 



where 

Fi(Ci) = InZi(Ci) = \nSpe^ Si = In 2 cosh 6 , 

Fr({r(})= In S P e^ = Fi° 2 T = 

= InZ^ ^ = F R ( Rl (i, R2 ( 2 , r 3 (3, R4C4) = 

= F R (C 1 , ~( 2 ,~( 3 ,~U), (5.4a) 

H({s}) = -V R ({s}) + J2 R f CfSf, 

feR 

Q = Cf + E R<pf> H/0 = q - R f <pf = Cf ■ ( 5 - 4b ) 

R 

The explicit form of cluster partition function is the following 
l -Z R = lz[%t = 2acosh(C~! - ~( 3 ) cosh(C~ 2 - C~ 4 ) + 

+dcosh(C 1 - C2 + Cs - L) + cosh(^ + ~C 2 + ~( 3 + C 4 ) + 
+26[cosh(^ + ~( 3 ) cosh(C 2 - C 4 ) + (5.5) 

+ cosh(C 1 - C 3 ) cosh(C 2 + C 4 )] , 

o = e"^; b = e~< iW - d = . 

Let us note that we neglect general factor exp{(Wi + AW + 2e)/8} in (5.5). 

It gives independent of temperature term in to free energy — /3 _1 ^ r ({C})- 
Let us introduce the following notations for derivatives 



,] (Ci) " m Ha f ' R ~^y 1234 ' 



CLUSTER APPROACH FOR ISING-TYPE MODELS 



85 



F?P,n = T^ g; /,/' = (1,2,3,4) = (R) . 

<■> [(f) <> {</>) 

The system of equations for averages rm = (Si) and fields Rtpi is obtained 
as in the previous section on the basis (2.5) and stationarity condition of 
T({C$) with respect to Ripi ( L mi = mi) 

m . = kjrW = FW(Ci) = tanhO , (5.7a) 
^i (1) (Ci) = 41; r = 0,l. (5.7b) 



Here index r numbers two clusters r = 0, 1 which contain bond 1 (see 
figure 2). Therefore there are two fields per one site. Then the system 
(5.7b) gives 8N equations in 8N unknowns Ripi where N is the number of 
A (or B) tetrahedron in the crystal. We have one equation (5.7b) for the 
field ip in the case of uniform field. The system (5.7b) is equivalent to the 
following system of equations 

pi(Si) = Sps 2r s ar s 4r PR(Si, S 2r , S 3r , S^) ■ (5.8) 

In particular from (5.8) the relation between intracluster CF follows: 

(S l 1 ) c pi =Fl l \c 1 ) = F«} t = (S[y pR . (5.9) 
One can obtain the relation for moment function from (5.8) 

M(R) = Sp l: 2,3,4PiP2P3P4e UR = Sp 2 ,3AP2p3P4e Uli (5.W) 

One can carry out the investigation of diagrams for F- function and = 

j^F in the higher order of cluster expansion similar as it was done in 

the case of two-particle cluster. This investigation leads to conclusion that 
FPCA gives the exact results for all characteristics of the system with the 
Hamiltonian (5.1) in the case of tree-like lattice. This lattice is constructed 
on the basis of tetrahedra. It can be formed from two-particle tree with 
z = 4 if all sites will change by tetrahedra. Now let us write some results 
in the homogeneous field. The free energy per tetrahedron in FPCA taking 
into account long-range interaction in MFA reads (m a = m): 

f = -^F{h) = '^^m 2 + 4kT ln2 cosh (3(-2kT In Z R (l), 
^Z R (C) = 2a + d + 46cosh2/?C + cosh4/3C, (5.11) 
The eigenvalue ^i(go) of the matrix J aa (qo) reads: 

Mqo) = J nWo) + Ji3(%) + 2Ji 2 (go) , (5.12) 
where % is the vector belong to the wave vector's star {q } (e — ► 0) 

{qo} = ±%x,±% y ; 9ox = (e,0,0); q 0y = (0, e, 0) . (5.13) 
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Fields (, ( contain the eigenvalue ^i(go) via molecular field v\(qo)m: 

( = (-2<p; ( = (-ip; ( = h + i/i(gb)m . (5.14) 
The system of equations for m, ip which follows from(5.7) reads 

m = tanh/3C = 4sinh2/?C [b + cosh2/3C] ^(C) • (5-15) 

Evidently from the first equation of (5.15)we can obtain the expression for 
ip = ip(m) and we can reduce (5.15) to one equation for parameter m. 

Further we shall use also the following intracluster CF (we present them 
only at Q = C) 

F W =m = tanh 0(; = 1 - m 2 , 

p t}\ - _ ,, -<* + cosh 4/?C 2 t>t c 1 -d 
^(0-^-2——^ m - 1 + 2a + 46 + d , 

P ^ _ _ +^-2a + cosh4/?C 2 t>t c 1 - 2a + rf 
F 2 (Q-F^ R -2 — ^ m - 1 + 2a + 46 + d > 

F l = ^fl = F 13,fl = ^il.i? = f uIr ' ( 5 - 16 ) 

5.2 Pair correlation functions of the reference system. 

The procedure for obtaining CFs of the reference system is similar to men- 
tioned above in the case of two-particle cluster. From (2.5) and (5.7) one 
has 

(S^) = Ff = if ^ , (5.17a) 

7(1) 5 ? x V" (!) (!) S 



Q = = Sij ~ Ri <p$ =(g ] + RVij • (5.17b) 



For obtaining equation for £jp we differentiate (5.7b) with respect to Q. 
This leads to relation (the notations R, i? r , 1, 2 r , 3 r , 4 r correspond to figure 
2) 

p( 2 ) JL p( n ) J 1 ) i p( U ) ,„(!) i p( U ) ,„(!) - 

*1 flVlj +*'l2,RR ( P2j + *13,RR<P3j + jf 'l4, J R^4j ~ 

- p (ll) p(H) p(H) fricA 

After cyclic permutation of the indices 1, 2, 3, 4, in (5.18) we shall find 
another three relations of (5.18) type. One can write for cluster Ri four 
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equation of the same type. Let us present these equations for cluster R r 
(r = 0, 1) in the matrix form 



Rr^lr),] = {Rr^if}', C(r),j = {Ci^} • 



7,(11) 



(5.19) 



(5.20) 



Here F^ 2 \ are the square matrices. They are formed by indices /, /' = 

l r ,2 r ,3 r ,4 r . The columnar matrices j^^ 1 ', £■ are formed by index j = 
1,2,3,4. 

In the following we consider only uniform case. Then on the basis of 
two matrix equations (5.19) (writen at r = and r = 1) and after some 
transformations we obtain the matrix equation 



/ 2cT } 



-A 



;(i) , ;(i) 

f(l) , 7(D 
**3j ^ *»3ij 

\Qj +C 4lJ ,/ 



(5.21) 



where 



/ F 1 F 2 F 1 
F^ = {F^S ff };F^ = g * g 

\Fi F 2 Fi 



It is necessary to obtain the closed equation for (ff from (5.21). Therefore 
we find the first element at the left hand side of (5.21) 



(5.22a) 
(5.22b) 



E 

reR,Ri 



A n = ^l 



11 



(5.23a) 



(5.23b) 



r-l 



Let us introduce matrices M 1 and M \ They are connected one with 
other as A and A^ 1 

M = (1 - A)" 1 = 1 + ( J p(2))-i j p(n). A = 1 - M- 1 , 

[M l )i r = {M~\ ri = (M- 1 )^; (M _1 ) = (IT 1 ),, . (5.24) 

Then the equation (5.23) for C^ 1 takes the form 

]T {<5 lr [l " 2(M _1 ) ] - (M _1 )i r } cg } = • (5.25) 

r£R 1 R 2 
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Taking into account (5.17) we write the equation for the reference CF P^ 2 ) 

-5ij . (5.26) 



E {Sir[-f^ j -2Po]-P lr }F 1 



(2) 
rj 



Here the matrix F is defined as 

p = M- 1 {F^)- 1 P- 1 = P (2) + P (11) . 



(5.27) 



Let us return to complex notations of bonds i = (nj,aj) for relation (5.26). 
We shall express the basic vectors of the Bravais lattice cell ai, a2, a3 in 
terms of vectors of the West's cell e a 

1 r i 1 r 

ai = - [-ei + e 2 + e 3 j ; a 2 = - [e x + e 2 + e 3 j ; a 3 = e 2 , 

bi = 2tt [-ei + e 3 ] ; b 2 = 2ir [ei + e 3 ] ; b 3 = 2ir [e 2 - e 3 ] . (5.28) 

Here b a are the basic vectors of reciprocal lattice. Then after Fourier trans- 
formation over n the relation (5.26) takes the form 



Ql 



1 



2P ]-P QQ1 (q)}^ ) Q ,(q; 



{(l-m 2 )- 1 -P(q)}F( 2 )(q) = -7. 
Here the matrix P(q) has the form 



P(q) 



/ 2P P 3 (q) Px(q) P 2 (q)\ 

P 3 (q) 2P P 5 (q) P 6 (q) 

Pl(q) P 5 (q) 2P P 4 (q) 

VP 2 (q) P 6 (q) P 4 (q) 2P / 



(5.29) 



(5.30) 



Pi(q) = 2P 2 cos y ; P 2 (q) = 2Pi cos -(q x - q y + q z ) , 

P 3 (q) = 2Pi cos ^(q x + q y + q z ); P 4 (q) = 2P 1 cos ^(q x + q y - q z ) , 



P 5 (q) = 2Pi cos -{q x - q y - q z ); P 6 (q) = 2P 2 cos . 



(5.31) 



The quantities Pq, P±, P 2 are expressed in the terms of functions F^ 2 \ Pi, 
P 2 with the aid of the relation 



Ip (q = ) = P= [p( 2 ) + p( n )] 



-i 



(5.32) 



Here q x , q y , q z are projections of dimensionless wave- vector q on the vectors 
e a and in the first Brillouin zone q a £ [— ir, tt]. From (5.32) we find 



P = V^U ; Pi = V~ 1 U l ;P 2 = V^U 2 , 
P + 2Pi + P 2 = (F^ + 2Pi + F 2 )~ l = - 



(5.33a) 



Zr 



Po - P 2 = (P (2) - F2)- 1 



1 Zr 



8 bc 2 + a ' 



8 C4 + 6c 2 — 0.5ZRm 2 ' 

(5.33b) 
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32 

V = —(bc 2 + a)(bc 2 + d)(c4 + bc 2 - 0.5Z r m 2 ) , 
Zr 

U = (36c 2 + a + 2ab + 2d) + (46 + 3d)6c 2 + ad - 

-0.5Z jR m 2 (2c 4 + 2bc 2 + d) - \z\ra 2 , 

8 

E7i = -(c4 - d)(6c 2 + a) + 0.5Z R m 2 (6c 2 + a) , (5.33c) 
U2 = (— &c 2 + a — 2d)c 4 + (2a — d)bc2 + ad — 

-0.5Z R m 2 (-bc 2 + a - 2d) - ^rTO 4 , 

8 

c 4 = cosh 4/3C; c 2 = cosh 2/3£ • 

Let us note that relations of (5.29) at T > T c were obtained firstly in [34]. 
The method developed in this work is acceptable only in the case T > T c . 
Moreover in [34] the Sleter model corresponding to the Hamiltonian (5.1) 
at $ = 0, was considered. For the mentioned above cases the matrix (5.30) 
coincides with corresponding matrix of [34]. 

The components of susceptibility tenzor (see next subsection ) will con- 
tain only certain eigenvalues of matrix P(q). Therefore we shall consider 
the procedure of diagonalizing of equation (5.29). Let us introduce nota- 
tion for eigenfunctions U /^(q) and eigenvalues e^(q) of matrix -P(q) — 2PqI. 
After unitary transformation we obtain 

- (2) „ , „ „ , , „ - (2) 

T (q) = ^+(q)^ 2 )(q)[/(q); ^ 2 )(q) = U(q)F (q)C/+(q) , 

^ 2 / ) ,(q) = E^(q)^ 2) (q)^/v(q) ) 

^ 2 )(q) = [e M (q) + 2P - (1 - m 2 )' 1 ]' 1 . (5.34) 
Here quantities e (U (q) are to be found from the equation 

e 4 - [£ P 2 (q)]e 2 - 2 [P^q) (P 2 (q)P 4 (q) + P 3 (q)P 5 (q)) + 
/=l 

+P 6 (q) (P 2 (q)P 3 (q) + P 4 (q)P 5 (q))] e + i?(q)i£(q) + 
+P 2 (q)P 2 (q) +P 2 (q)P 2 (q) - 2P 2 (q)P 3 (q)P 4 (q)P 5 (q) - 
-2Pi(q)P 6 (q) [P 2 (q)P 5 (q) +P 3 (q)P 4 (q)] =0. (5.35) 

We consider such values q for which the problem of eigenfunctions and 
eigenvalues due to the symmetry of -P(q) may lead to solution the equations 
of second order. For the values q = (q, q, q z ) we find 



2ei, 2 (q) = P 3 (q) +P 4 (q) ± Jv+(q). 



2e 3 ,4(q) = -P 3 (q) - Pi(q) ± V^Co) > ( 5 - 36a ) 
P±(q) = [P 3 (q) - P 4 (q)] 2 + 4 [P x (q) ± P 2 (q)] 2 ; 
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C/(q) 



( Wi(q) ^(q) ^ n (q) 

Wi(q) ^12 (q) -^n(q) 

¥>i2 (q) -vii(q) ^(q) 

V^ 2 (q) -vn(q) -V12 (q) 



^(q) \ 
-^(q) 
-^Ti(q) 
^n(q) y 



(5.36b) 



vii(q) = -¥'22(q) = -s[i + *^(q)] 



1/2 



vi2 (q) = 



V2 

[l-Srsid)] 1 / 2 



r 5 (q 
For the values q 



V2 
1 



(5.36c) 



1/2 



±P|(q) ;fl 4 (q) 



P 3 (q) ~ P 4 (q) 
Pl(q)+5P 2 (q) 



(?x> 9y) 0) we obtain 



2ei, 3 (q) = Pi(q) +P 6 (q) ± >Jv+(q), 
2e 2 , 4 (q) = -Pi(q) - P 6 (q) ± V^fa) , 



P 3 ± 2 (q) = [P 1 (q) 

/Wi(q) 
^(q) 
^n(q) 
V^i2(q) 



f/(q) 



vii(q) = -^22 (q) = 



- P 6 (q)] 2 +4[P 3 (q) ± P 2 (q)]- 

W2(q) ^n(q) ^(q) \ 

^n(q) ^i2_(q) 

Vi2(q) -^n(q) 

-^n(q) -<Pi2 (q) 

[l + tfr^q)] 1 / 2 



-^n(q) 
-^12 (q 
^n(q) J 



(5.37a) 



(5.37b) 



¥>i2 (q) 



^<j(q) 



[l-Srsfr)] 1 ' 2 



V2 



V2 



(5.37c) 



P 5 (q) 



i + 4^(q) 



1/2 



; P 5 (q) 



Pi(q)-P 6 (q) 
P 3 (q)+5P 2 (q) 



Let us note that in [34] the eigenvalues were obtained only for vector 
q = (q, Qi Qz) and the long-range interaction was not considered. 

5.3 Susceptibility of KD 2 PO i system. 

The tensor of isotermical susceptibility of KD2PO4 system can be expressed 

T (2) 

in terms of the pair CF T a p (q) taking into consideration long-range inter- 
action in MFA in the following way 



Xafe(q) = d *adf3b L FaP (<l) 

-1 



^ 2 )(q)= ^(q 



pj(<c 



d az — d c ] 
d y = d2 y - 



dx — di x — —ds x ; d2x 
-d4 y ; di y = —dz y = . 



d^x = , 



(5.38a) 



(5.38b) 
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Here d aa is an effective dipolar moment of hydrogen bond. It contains 
the contribution of ionic subsystem via spin-phonon interaction. We write 
the longitiudinal and transversal components of tenzor Xab(q) hi terms of 

~ (2) T (2) 

eigenfunctions W a/X (q) and eigenvalues T\i (q) of the matrix T a p (q) • 

Xx M) = 4 E^i 2) (q) - w zMf , (5-39) 

X W (q) = 4E^(q) [^(q) - ^(q)] 2 . 

The symmetry of matrix ^^ 2 ^(q) which is related with the symmetry of 
matrix P(q) at arbitrary value of q does not coincide with the symmetry 
of matrix J a p(q). Here we consider only such values of q when P(q) and 
/(q) has equal eigenfunctions. Then on the basis (5.38a) and (5.34) one can 
write 

4 2) (q) = ["(I " m 2 )- 1 + 2P + e M (q) - ^(q)^' 

^(q) = E ^(q)^«'(q)^>(q) - 2P , 



^(q) = E H/ "^(q) J "«'(q)w / "a>(q)- 



(5.40) 



a,a' 



We present here several cases of application of relation (5.40). At q — > we 
use the expressions for J Ql g(q) from [59] 
1- q = q z = (0,0, q z ) 



J{q z ) = 



( J\i{q z ) J\2{q z ) J\z{q z ) J\2{q z )\ 

Ju(q z ) J\i{q z ) Ji2(q z ) Ji2,{q z ) 

Jis(q z ) J\2{q z ) Jn{q z ) J\2{q z ) 

\Ji2(q z ) Jn(q z ) Jn{q z ) Jn(q z ) J 



Ju(qz) >° Ju + a z ; 
Jis(q z ) q ^~^ Jis — a z , 



W(q z ) = U(q z ) 



Ju(q z ) q ^-^ J12 , 




vi,3(<lz) = Jn{q z ) + Ji3(q z ) ± 2Ji 2 (g z ) Q£ -^ Ju + J13 ± 2J i2 , 



"2,4(q z ) = Ju(q z ) ~ Jiz{q z ) q£ ~> Ju ~ J13 ± 2a z , 
ei >2 (fe) = 2[P 2 ± P 1 cos ^]; e 3A (q z ) = -2P 2 , 

(24) 2 

Xaa(g^) 

2(2d a ) 2 



(5.41a) 

(5.41b) 
(5.41c) 

(5.41d) 



= {2(P + 2Pi + P 2 ) - 8Pi sin 2 f - (1 - m 2 )^ 1 - ^i(fc)}- 1 , 
= {2(P - P 2 ) - (1 - m 2 )- 1 - v 2 {q z )Y 1 ; a = x,y . (5.41e) 
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2. q = q= (e,e,0), e^O 



m = 

Ju(q) 



(JiM 

Jl2(q) 

\Ju{q) 
Jll ~\~ dxy 'i 

J12 - a 



xy-i 



*-xy i 



W(q) = U(q) 



Jii{q) J i3(q) J u(q)\ 

Jn(g) Ju(q) Ji3{q) 

Ju{q) Ju(<i) Ji 2 (g) 
Jis(q) J\2{q) Ju(q)J 

J\2{q) = J\2 + O-xy , 

Ji3(q) = Ji3 - a 

1 1 1 

1 -1 -1 

-1 1 -1 

-1 -1 1 

vi,M = Ju(q) + ^13(9) ± [^12(9) + Ju(q)\ 
^2,4(9) = J11 (?) - ^13(9) ± [^12(9) - ^14(9)] , 
^2(9) = -hi - J13 + 4a XJ/ ; 1/4(9) = Jn ~ Ji3 
ei )3 (g) = ±4Pi + 2P 2 ; €2,4(9) = "2P 2 , 

Xzzjq) 

Xaa(q) 



(5.42a) 

(5.42b) 
(5.42c) 



Jll + Jl3 ± 2Ji2 , 



(5.42d) 



(24) 2 



{2(P + 2Pi + P 2 ) - (1 - m 2 )" 1 - ^(g)}- 1 

- ^r 1 • 



x: {2(p - p 2 ) - (1 

M=2,4 



m 



(5.42e) 



3. q = o, x= (e,0,0), y = (0,e,0), e^O 



J(a) 



/ Jn (a) Ji 2 (a) Ji 3 (a) Ji 2 (a)\ 

Ji 2 (a) J 22 (a) Ji 2 (a) J 2 4(a) 

Ji3(a) Ji 2 (a) Jn (a) Ji 2 (a) ' 

\Ji 2 (a) J24(a) Ji 2 (a) J 22 (a)/ 

Jn (a) = Jn + 6(a); J 22 (a) = J11 + a(a) , 
Ji3(a) = J13 - b(a); J 2 4(a) = J13 - a(a) , 
J 12 (a) = J i2 ; 6(a) = a (a) , 



VF(a) = 

^1,3(9) 

^2(9) = 
1/4 (g) = 

£1,2(9) = 
Xzzjq) 
2(24) 2 

Xaa(g) 

2(24) 2 



Ilia) = \ 



1 -1 
1 1 

M -1 



V2 


-V2 




\ 






-^2/ 



(5.43a) 



(5.43b) 



(5.43c) 



= Jn (a) + Ji3(a) ± Ji 2 (a) = J n + J i3 ± 2J i2 
Jn (a) - Ji 3 (a) = J11 - J13 + 26(a) , 
J 22 (a) - J24(a) = Jn - Ji 3 + 2a(a) , 
= ±4P! + 2P2; e 3 , 4 (9) = "2P 2 , 

{2(P + 2Pi + P 2 ) - (1 - m 2 )" 1 - Vl (q)Y 



(5.43d) 



(5.43e) 



{2(P - P 2 ) - (1 



m 



^(9)}" 



A* 



2, a = x 
4, a = y 



More detailed investigation of ^"-dependence for susceptibility of KD2PO^ 
system will be carried out in the separate paper. 
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6 Conclusion 

Thus in the present work we have suggested the cluster method for compu- 
tation of the correlation functions of the Ising-type systems with short-range 
interactions and with arbitrary value of spin. The method is based on the 
calculation of generating function (that is the logarithm of partition func- 
tion for a system in nonuniform external field) . Here we restricted ourselves 
only to the first order of cluster approximation for generating function. The 
investigation of some models (S z = ±1, S z = —2,0, +2) has been carried 
out within two-particle cluster approximation. The static susceptibility of 
KD2PO4 crystal has been obtained within four-particle cluster approxima- 
tion. 

The obtained here form of equations for CFs is similar to the form which 
one can obtain within random phase approximation (molecular field approx- 
imation for generating function) . Moreover in the case of Bravais lattice the 
q-dependence of CFs is similar to one obtained within random-phase ap- 
proximation. Nevertheless within cluster approximation we obtain the more 
precise expressions for the temperature depending coefficients. In the case 
of tree-like lattices, in particular, for ID lattices, the cluster approxima- 
tion gives the exact expression for all static characteristics of the Ising-type 
systems. The more detailed study including numerical calculations of ex- 
pressions obtained in the present work will be given elsewhere. 
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